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Microwave holography has proven to be a powerful technique for various evalua- 
tions, diagnostics , and RF performance improvements for large reflector antennas. 
This technique utilizes the Fourier transform relation between the complex far-field 
radiation pattern of an antenna and the complex aperture field distribution. Re- 
sulting aperture phase and amplitude distribution data can be used to precisely 
characterize various crucial performance parameters , inciuding panel alignment , 
panel shaping, subreflector position , antenna aperture illumination , directivity at 
various frequencies , and gravity deformation effects. 

The methodology of the data processing presented in this article was developed at 
JPL and has been successfully applied to the DSN 34-m beam waveguide antennas. 
A companion article in this issue (Rochblatt and Seidel) describes the application 
of this technology to the DSS-13 antenna. The antenna performance was improved 
at all operating frequencies (wide-bandwidth improvement) by reducing the main 
reflector mechanical surface rms error to 0.43 mm. At Ka-band (32-GHz), the esti- 
mated improvement is 4.1 dB, resulting in an aperture efficiency of 52 percent. The 
performance improvement was verified by efficiency measurements and additional 
holographic measurements. 


I. Introduction 

Microwave holography is a measurement technique that 
has now been successfully applied to improve the per- 
formance of many types of reflector and array antennas 
throughout the world [1-7]. 1 The raw data (the observ- 
able) for this technique is the complex far-field pattern 

1 B. L. Seidel and D. J. Rochblatt, DSS-13 Beam Waveguide An- 
tenna Phase I Final Report, JPL D-8451 (internal document), Jet 
Propulsion Laboratory, Pasadena, California, pp. 4.1-4.23, May 15, 
1991. 


of the antenna under test. For large reflector antennas, 
geostationary satellites are commonly used as convenient 
far-field signal sources. The beacon continuous wave (CW) 
signal from a geostationary satellite is usually used with 
a narrow-bandwidth receiver, while a transponder signal 
is usually the preferred choice for a wide-bandwidth re- 
ceiver. Figure 1 presents a typical narrow-bandwidth re- 
ceiver system architecture based on an HP 8510B analyzer 
and an external phase-locked loop (PLL) (with a variable 
PLL bandwidth of 1 to 20 Hz). The reference antenna is 
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needed as a phase reference for the measurement, as well 
as to keep the narrow-bandwidth receiver in phase lock 
with the carrier [8,9]. 

This article provides a detailed analysis of the unique 
methodology and data processing developed at JPL by the 
author. 


1/2APL = l/2[P f P -b PQ] 
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Phase(APL) = —e cos (j> 


II. Mathematical Algorithms an d 

The mathematical relationship between an antenna 
far-field radiation pattern (T) and the antenna surface- 
induced current distribution (J) is given by the exact ra- 
diation integral [10] 



( 5 ) 


T (u, v ) = J J J (x / ,y')exp ; ^'[exp--’* i '( 1 - cos9 )] 
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and for a paraboloid: 



( 6 ) 


where 

z(x f , y f ) = the surface s 

(u, v) = directional cosine space 
6 = observation angle 

For small observation angles (9 = 0 deg) of the far- 
field pattern, this expression becomes a Fourier transform 
relationship 

T (u y v) = J J J (x f , y')exp jkz exp 7 *^* +,jy ^ dx'dy* 

8 

( 2 ) 

To derive the residual surface error, geometrical optics- 
ray tracing is used to relate the small normal error, e, 
directly to an aperture phase error in a main reflector 
paraboloid geometry (Fig. 2). In this small error approx- 
imation, it is assumed that the aperture phase error is 
entirely due to the projected surface errors. In addition, 
the three-dimensional structure of the surface error is not 
recovered; rather, an axial component equal to the aver- 
age error over the resolution cell size is recovered. The 
normal error can then be computed using the paraboloid 
geometry. In practical application, when many resolution 
cells resolve a single panel, high accuracy (typically A/250 
at 12 GHz) can be achieved in determining panel setting 
screw adjustments. It can be computed from the geometry 
in Fig. 2 


Allowing for the removal of a constant phase term and 
substituting Eq. (4) into Eq. (2) yields: 

T(u,v) =exp ~i 2kF J J [\J(x\y , )\exp^ 7r * cos< ^ exp jkz '] 

s 

x expi k ( ux +vy ) dx'dy' (7) 

where [...] = J e (equivalent current distribution). 

For the processing of sampled data, the associated dis- 
crete Fourier transform (DFT) is utilized: 

Nl/2-1 N2f2-l 

T(pAu } qAv) =sxsy ^ ^ J(nsx,msy) 

n = — Nl /2 m = — 7V2/2 

x exp J ’ 2,r ($T + # 2 ) (8) 

where 

N 1, N 2 = the measured data array size 

sx,sy= sampling intervals in the aperture coor- 
dinates 

n, m, p, q = integers indexing the discrete samples 
Au y Av = sampling intervals in u,v far-field space 
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Since the magnitude of the far-field pattern is essen- 
tially bounded, the fast Fourier transform (FFT) is usu- 
ally used for computation, and is symbolized here by (F). 
Solving for the residual normal surface error and substi- 
tuting Eq. (6) yields 

c ( x ’ y ) = ^\A + ^4^- phase {exp-' 2 **’ F^pXu.w)]} 

( 9 ) 

The spatial resolution in the final holographic maps is 
defined at the 50-percent level of the Fourier transform 
of the truncation function [7,11]. Rectangular and trian- 
gular windowing are commonly used, with competing ad- 
vantages to both. For rectangular windowing, the lateral 
resolution is derived [7,11] 


where 

D = main reflector diameter 

N — the square root of the total number of data points 
k — sampling factor, usually 0.5 < k < 1.0 

The lateral resolution is inversely proportional to the num- 
ber of sidelobes measured. For the DSN 70-m antennas, 
the required resolution of 0.4 m is achieved with a data 
array size of 197 x 197. For the DSN 34-m antennas, the 
required resolution of 0.32 m is achieved with a data array 
size of 127 x 127 [12]. 

Figure 3 shows the complex Fourier transform relation- 
ship between the far-field and the aperture functions for 
a 34-m antenna (DSS 13). The far-field amplitude and 
phase, Figs. 3(a) and 3(b), respectively, are measured on 
rectangular coordinates of 127 x 127 with sampling in- 
tervals of 34.8 mdeg (the sampling factor is 0.84). Fig- 
ures 3(c) and 3(d) show the aperture amplitude and sur- 
face error function with a lateral resolution of 0.32 m. 

The accuracy in each resolution cell of the final holo- 
graphic maps was formulated from a simulation study [9], 
which agrees well with analytically derived expressions [2]. 
In this study [9], Eq. (11) is derived 


a = standard deviation (accuracy) in recovering the 
mean position of a resolution cell 

A = wavelength 

SNR = beam peak voltage signal-to-noise ratio in the 
test (antenna) channel 

The accuracy across the holographic map varies with 
the aperture amplitude illumination. The accuracy is bet- 
ter at the center of the dish and gradually becomes worse 
toward the edge, where the illumination falls off rapidly [9]. 
For a uniformly illuminated antenna design, the accuracy 
remains relatively constant, becoming quickly worse just 
at the edge, where the illumination falls off rapidly. Typi- 
cally, accuracies of 0.05 mm to 0.10 mm were achieved with 
a corresponding SNR of 73 dB to 60 dB. Equation (11) 
also indicates that the accuracy is inversely proportional 
to the resolution. This is due to the larger averaging area 
available at the larger resolution cell, as expected. 

The resulting aperture phase function (the dirty map) 
needs to be corrected for modulo-27r phase errors. These 
occur due to the small measurement wavelength and a 
large phase error that is partially due to pointing and 
subreflector position errors. Figures 4 and 5 are exam- 
ples of aperture phase maps before and after correction of 
modulo-27r phase errors, respectively. 

The next phase cleanup algorithm removes from the 
aperture function phase error components that are due to 
pointing and subreflector position errors. This can be done 
by fitting polynomials to the aperture phase data whose 
orders are equal to the expected order of the errors. Either 
the Zernike [13] or the modified Jacobi polynomials [14] 
can be used. Determination of pointing error-linear term, 
astigmatism, focus shift, coma, and spherical aberration 
are readily identified. The orthogonality properties of the 
Jacobi polynomial allow solving for its coefficients by 

1 2jt 
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where 


e(s l ,$) = surface error function 


( n = Neumann number 
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In addition, the above polynomials can be used to sepa- 
rate systematic (low-frequency) errors from random (high- 
frequency) error components by filtering (Fig. 6) [6]. 

Another technique for removing phase errors due to 
pointing and subreflector position errors has been devel- 
oped via the global best fit paraboloid. Here, the entire 
data set is weighted-least-squares fitted to the paraboloid, 
permitting six degrees of freedom in the model, three ver- 
tex translations, two rotations, and a focal length change. 

Minimization is performed on the sum squares of the 
residual path length changes 5, given by 

TV 3 

S = £r(APL < ) 2 ^ (14) 

1=1 


where 

T = support domain masking operator 
A PLi — path length change 

Ai = amplitude weighting factor 

The minimum for 5 is found by solving the six partial 
differential equations simultaneously 
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X 0 ] 

Yo > = vertex coordinates translation 

J 

a = rotation about the x-axis 
) 3 = rotation about the y-axis 


and 



where 

F — focal length of the original paraboloid 
F 9 = focal length of the best-fit paraboloid 


The antenna surface axial rms error is then computed with 
respect to the position of the fitting paraboloid 


RMS - 



r(i/2APL,) 2 



(17) 


The subreflector position error is found via the six pa- 
rameters computed above. Figure 7 presents a large linear 
phase error term due to pointing errors where errors due 
to subreflector position are minimal. 

It is correct to apply the best-fit paraboloid algo- 
rithm to either the conventional Cassegrain paraboloid- 
hyperboloid or dual-shaped reflector systems, even though 
the latter do not use a paraboloid as the main reflector. 
Either design is a plane-wave-to-point source transformer, 
differing only in the field intensity distribution. 

The resultant aperture function at the end of this pro- 
cess is defined here as the effective map since it includes all 
phase effects that are contributing to the antenna perfor- 
mance [15]. These frequency-dependent effects include the 
subreflector scattered (frequency-dependent) feed phase 
function and struts’ diffraction effects. Removal of the feed 
phase function and subreflector support structure diffrac- 
tion effects results in a frequency independent map, which 
is defined here as the mechanical map. 

By deriving panel adjustment based on the effective 
map, the surface shape will conjugate to the phase errors, 
optimizing the performance of the antenna at a single fre- 
quency equal to the measurement frequency, while degrad- 
ing the performance at all other frequencies. For antennas 
operating at a single frequency, this procedure is advan- 
tageous. However, many antennas (including the DSN’s) 
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operate at several different frequencies and require wide- 
band width performance. For these antennas, the mechan- 
ical map must be used to derive panel setting information. 

Figure 8 is a feed phase function derived for the DSS-13 
antenna, a 34-m beam-waveguide antenna at 12.198 GHz, 
using Geometrical Theory of Diffraction (GTD). Figure 9 
is an effective map derived for the same antenna at 
12.198 GHz. Subtracting the feed phase function from 
the effective map and masking the struts’ diffraction ef- 
fects resulted in the mechanical map presented in Fig. 10. 
Successful removal of the frequency-dependent feed phase 
function must result in reduction of the rms error for the 
mechanical map. This is due to the uniqueness of such a 
solution. As shown above, the rms decreased by 27.5 per- 
cent, from 0.98 mm to 0.71 mm. 


Panel setting information is derived by sorting all the 
data points within each panel and performing a rigid-body 
least-squares fit. The algorithm allows for one translation 
and two rotations (S k , a k , and f3 k ), hence, a rigid body 
motion (Fig. 11). For each panel and its associated n data 
points, the motion parameters are solved via Eq. (18) 
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This mathematical process increases the accuracy in de- 
termining the screw adjustment correction by a factor of 
yfn, since all n data points within each panel are used to 
derive each screw adjustment. The three parameters from 
the least squares fit above are also used to recompute the 
predicted surface errors within each panel that will result 
after panel setting. Figures 12, 13, and 14 are, respec- 
tively, the before, predicted, and after mechanical maps 
derived during holographic measurements of DSS 13 and 
applying the above algorithms. The antenna was found to 
have an rms surface error of 0.88 mm, which was reduced 
to 0.45 mm. Figure 13, the predicted surface map (rms = 
0.36 mm), indicates that panels in rings 8 and 9 were de- 
formed, which was then confirmed in Fig. 14, the derived 
surface after panel setting. 

The Appendix by S. Stewart, D. Rochblatt, and B. Sei- 
del of the companion article in this issue describes in 
detail the cause for these deformed features due to us- 
ing DSS- 15 main reflector panels on the DSS-13 beam- 
waveguide (BWG) antenna. The reduction in rms from 
0.88 mm to 0.45 mm was achieved within one session of 
panel setting, which limited the screw adjustments to the 
nearest 1/8 of a turn. 2 Figure 15 is a sample printout 
of the panel-setting screw-adjustment listing produced by 
this software package. 

Another significant feature of the mechanical map is 
that surface tolerance efficiency can be computed at fre- 
quencies other than the measured one. This is possible 
since the aperture phase function has been cleaned of all 
effects that are frequency dependent. Surface tolerance 
efficiency can be computed by integrating the entire aper- 
ture function. In this computation, it is assumed that 
the aperture amplitude function is also frequency inde- 
pendent. The introduced error in this assumption is thus 
very small. 
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where 


<f) mt — phase value at the measured frequency 
sub k = any other frequency 


An effective surface error e can be computed by equating 
the above expression to that of Ruze [11,16] 


2 ibid. 
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t] = exp — (47r£:/A) 2 (20) 

Figure 16 shows the estimated surface loss versus fre- 
quencies of DSS 13 before and after panel setting, indicat- 
ing a performance improvement of 4.1 dB at 32 GHz, which 
resulted in an aperture efficiency of 52 percent. Figure 17 
shows the measured efficiency of the DSS-13 antenna be- 
fore and after panel setting as a function of elevation angle 
at X-band (8.45 GHz). As shown in Fig. 17, the peak an- 
tenna gain increased by 0.21 dB and moved from the anti- 
symmetric 57-deg elevation to the desired 45-deg elevation 
that was specified as the required antenna rigging angle. 

Figure 18 is a gravity deformation of DSS 13 between 
the elevation angles of 12.7 deg and 46.5 deg. Gravity 


deformation maps are highly important in deriving accu- 
rate structural deformation models for large reflector an- 
tennas. It is observed from this image (Fig. 18), that the 
bypass shroud on the right side of the antenna is exert- 
ing a force in the outer axial direction, creating the ob- 
served signature. It was found that the disturbing shroud 
forces are attributed to eccentricity of the shroud axes of 
rotation with respect to the antenna elevation axis. 3 The 
eccentricity is also confirmed by dial-indicator measure- 
ments. Figure 19 presents the major software components 
of the microwave holography methodology in a block dia- 
gram. 


3 R. Levy, DSS-13 Antenna Structure Measurements and Evalua- 
tion , JPL D-8947 (internal document), Jet Propulsion Laboratory, 
Pasadena, California, October 1, 1991. 
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Fig. 3. Complex Fourier transform relating the antenna far-field and aperture functions. 









Fig. 4. Phase map before correction for modulo 27T. 
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Fig. 5. Phase map after correction for modulo 27r. 
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Fig. 6. Modified Jacobi polynomial decomposition separates structural deformation from panel setting errors. 






Fig. 7. Phase map due to pointing and subreflector position errors. 



Fig. 8. Feed phase function for DSS-13 34-m antenna at 12.198 GHz. 
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Fig. 9. Effective map. 



Fig. 10. Mechanical map. 
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Fig. 12. DSS-13 mechanical surface as found. 
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Fig. 13. DSS-13 computer-predicted mechanical surface based on methodology. 



Fig. 14. DSS-13 mechanical surface after panel setting. 
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SCREW ADJUSTMENTS FOR RING 1 


PANEL 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


A 

0.0464 

0.0296 

0.0172 

0.0473 

0.0296 

0.0507 

0.0593 

0.0808 

0.0374 

0.0654 

0.0306 

0.1012 


B 

0.0309 

0.0047 

0.0071 

0.0096 

-0.0044 

0.0213 

0.0232 

0.0837 

0.0089 

0.0300 

0.0230 

0.1022 


C 

0.0171 

-0.0006 

0.0146 

0.0084 

0.0179 

0.0110 

0,0303 

0.0441 

0.0229 

0.0258 

0.0261 

0.0494 


D 

0.0065 

-0.0047 

0.0204 

0.0075 

0.0350 

0.0031 

0.0358 

0.0135 

0.0337 

0.0226 

0.0284 

0.0087 


E 

0.0325 

0.0243 

0.0248 

0.0461 

0.0521 

0.0403 

0.0665 

0.0407 

0.0515 

0.0611 

0.0337 

0.0478 


F 


SCREW ADJUSTMENTS FOR RING 2 
BCD 


PANEL 

1 

2 

3 

4 

5 

6 


A 

0.0412 

0.0150 

-0.0050 

0.0009 

0.0114 

0.0188 


0.0224 

0.0242 

0.0027 

0.0012 

0.0116 

0.0071 


0.0077 

0.0313 

0.0086 

0.0015 

0.0117 

- 0.0020 


0.0142 

0.0020 

0.0038 

0.0099 

0.0138 

0.0087 


E 

0.0275 

0.0009 

- 0.0012 

0.0079 

0.0132 

0.0156 


F 

0.0446 

-0.0004 

-0.0075 

0.0054 

0.0125 

0.0244 


Fig. 15. Panel-setting screw-adjustments (inches) computer listing for DSS-13 34-m 
antenna derived from scan JPL106 127 X 127. 



Fig. 16. Estimated DSS-13 surface loss. 



EFFICIENCY 



AUG/SEPT 1990 


WITHOUT ATMOSPHERE 


ELEVATION ANGLE, deg 


Fig. 17. DSS-13 X-band (8.45-GHz) efficiency at the Cassegrain focus before (bottom) and 
after (top) application of the holography methodology. 


Fig. 18. Gravity deformation signature derived for DSS 13 from measurements at two 
elevation angles (12.7-46.5 deg). 
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SIMULATION CAPABILITY 


DATA-PROCESSING AND DISPLAY CAPABILITY 



Fig. 19. Simplified block diagram of the microwave holography software methodology. 
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